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jü 3E: 目前 , 通过 卫星 遂 感 、 低 空 适 感 影像 自动 ,半自动 化 解 译 梯田 信息 有 一 定 的 研究 进展 ， 但 是 受 数据 获取 成 本 、 
精度 、 解 译 方法 单一 等 限制 ， 只 限于 大 面积 提取 梯田 区 域 。 低 成 本 地 进行 梯田 的 田 面 精准 提取 以 及 面积 统计 仍 需 进 一 
步 研究 。 基于 面向 对 象 方法 分 别 对 0.5m 分 辨 率 的 无 人 机 正 射 影像 和 地 形 指数 ,及 两 种 数据 的 结合 进行 梯田 区 域 分 割 、 
提取 及 面积 统计 。 结 果 表 明 ， 将 正 射影 像 数 据 与 地 形 指 数 结合 的 梯田 的 田 面 提取 结果 优 于 基于 单一 数据 源 。 

关键 词 : 梯田 ; 无 人 机 正 射 影像 ; 地 形 指数 ; 面向 对 象 ; 监督 分 类 
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Extraction of terrace information based on UAV image and topographic index 
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Abstract: At present, semi-automatic interpreting terrace information had certain research progress through satellite remote 
sensing and low-altitude remote sensing imagery. However, limited by data acquisition cost, accuracy, and single interpretation 
methods, it was currently limited to large-scale extraction of terraced areas. Accurate extraction of terrace field and area statistics 
at low cost still require further study. In this study, based on the object-oriented approach, the orthophoto and terrain index of 
UAV with 0.5m resolution, and the fusion of the two kinds of data were used for extraction of terrace field and area statistic. 


The results show that with the fusion of orthophoto and terrain index it is superior to those based on a single data source. 


Key words: terrace; UAV orthophoto; topographic index; object-oriented ; supervised classification 


c 应 用 于 梯田 信息 提取 。 张 雨 果 品 , 党 恬 敏 四 等 研究 者 基于 分 辩 率 
一 相对 较 低 的 单一 卫星 遥感 数据 源 ， 试 错 法 找到 最 优 尺 度 〈 主 观 
梯田 的 修筑 是 人 类 在 坡 耕 地 上 探索 的 一 种 水 土 保 持 措 施 ， 性 强 , 缺乏 定量 分 析 ), 采用 了 基于 规则 ( 需 人 为 设 定 规 则 ， 主 
很 大 程度 上 开发 了 坡 耕 地 的 农业 生长 潜力 ， 具 有 鞭 水 、 保 土 作 观 性 强 ) 和 最 邻近 (K-nearestneighbor, KNN) 的 面向 对 象 分 
用 巾 。 长 期 以 来 ， 通 过 目 视 解 译 是 获取 梯田 信息 的 主要 技术 手 。 ”类 方法 。 在 分 类 之 后 手工 参与 错 分 ， 漏 分 调整 ， 属 于 半自动 化 
段 。 但 限于 经 验 知识 ， 技 术 成 本 等 ， 对 于 梯田 的 数量 、 面 积 梯田 区 域 提取 。Eckert 等 人 中 结合 4 波段 的 多 光谱 数据 和 DSM 
分 布 等 信息 的 调查 资料 不 完善 ， 导 致 梯田 在 水 土 保持 中 的 贡献 数据， 运用 面向 对 象 分 类 方法 梯田 进行 提取 ， 将 光谱 信息 和 
缺乏 定量 研究 。 随 着 高 分 辨 率 影像 获取 技术 和 处 理 技术 的 不 断 ” DSM 结合 相 较 于 单独 使 用 光谱 信息 提取 精度 更 高 。 与 卫星 数据 
发 展 ， 基 于 纹理 特征 的 傅 里 叶 变 换算 法 ， 模 板 匹 配 算 法 等 梯田 相 比 ， 无 人 机 航 摄 系统 有 低 成 本 ， 高 效率 快速 及 时 获取 高 分 辨 
言 息 提取 成 为 可 能 1。 但 由 于 纹理 特征 不 明显 ,不 规则 ,以 及 ”” 率 影像 的 优势 。Han 等 人 外 基于 DEM， 和 基于 DEM 与 正 射影 
傅 里 叶 变 换算 法 窗口 大 小 影响 ， 模 板 匹 配 算法 的 模板 尺寸 固定 。 像 与 红外 波段 结合 的 两 类 数据 源 ， 对 比分 析 基 于 规则 梯田 分 类 
不 灵活 等 限制 ， 大 量 漏 提 和 和 错 提现 象 。 而 且 算 法 只 考虑 梯田 的 结果 。 但 是 该 研究 并 未 分 析 单 独 基于 正 射 影像 与 红外 波段 结合 
纹理 特征 和 灰 度 特征 ， 对 于 梯田 的 颜色 ， 形 状 ， 坡 度 等 地 形 指 ”的 数据 源 的 分 类 情况 , 缺乏 一 定 的 严谨 性 和 说 服 力 。 Diaz-Varela 
数 并 未 充分 利用 。 SEDI, 结合 DEM 与 多 光谱 影像 数据 , 通过 面向 对 象 的 分 析 方 法 
面向 对 象 分 类 方法 能 充分 利用 梯田 的 颜色 ， 进行 梯田 区 域 的 提取 ， 但 研究 区 域 田 面 植被 覆盖 率 较 高 ， 提 取 
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录用 入 BAUR, $: 无 人 机 影 hips yA 
精度 受到 一 定 影响 。 筑 物 以 及 道路 ， 田 面 无 覆盖 物 ， 田 块 形 状 特征 较 短 ， 较 宽 ， 梯 

车 于 面向 对 象 分 析 ， 结 合 高 分 辩 率 的 影像 数据 和 DEM 数 1 数量 较 少 。 区 域 二 的 梯田 区 有 部 分 建筑 物 ， 田 面 有 部 分 庄稼 
据 ， 对 大 面积 梯田 区 域 提取 研究 已 有 一 定 进展 。 但 是 存在 如 下 积 物 ， 田 块 形 状 特征 整体 舒 长 ， 较 窗 ， 田 块 边缘 曲线 光滑 ， 
问题 : a) 几乎 所 有 研究 者 都 停留 在 大 面积 提取 梯田 区 域 ， 田 面 | 块 数量 较 多 。 区 域 三 的 梯田 区 有 部 分 建筑 物 ， 田 面 履 有 大 面 
精准 提取 ， 面 积 统计 等 信息 研究 目前 几乎 没有 研究 进展 ， 没 有 积 积 雪 ， 田 块 形状 蚁 归 曲 折 。 三 块 研究 区 基本 代表 了 样 区 的 梯 
满足 精准 农业 的 需求 ,而 且 大 部 分 数据 源 为 卫星 遥感 影像 ， 获 取 1 地 表 特 征 。 


难度 大 ， 成 本 昂贵 ;b〉 目 视 解 译 获 得 多 尺度 分 割 的 最 佳 分 割 尺 
E, 停留 在 定性 阶段 , 缺乏 定量 分 析 ;c) 分 类 方法 单一 ， 且 对 于 
精度 评价 几乎 都 未 提 及 评价 样本 的 获取 。 评 价 样本 越 客 观 ， 评 
寺 果 越 具 有 参考 性 。 

针对 目前 梯田 提取 存在 的 问题 ， 本 文 结合 无 人 机 影像 低 成 
本 ， 高 效率 快速 及 时 获取 的 优势 ， 应 用 无 人 机 拍摄 的 甘肃 龙泉 
0.5 m 高 分 状 率 正 射影 像 影像 数据 ， 对 无 植被 覆盖 的 水 平 梯田 
进行 提取 方法 研究 。 采 用 面向 对 象 分 类 技术 ， 对 大 面积 梯田 区 
域 的 田 面 进行 提取 ， 并 统计 相应 梯田 面积 。 首 先 ， 对 无 人 机 拍 
摄 的 正 射 影像 数据 进行 预 处 理 和 DEM 数据 地 形 指数 分 析 ;， 其 


AS 


次 应 用 胡 忠 文 等 人 Uo 区 的 最 优 尺 度 研 究 方法 ， 基 于 区 域 合 并 

演化 分 析 选 择 评 估 最 佳 分 割 尺度 ， 然 后 进行 分 类 特征 筛选 ， 最 图 1 甘肃 省 榆 中 县 龙泉 乡 DEM 和 正 射影 像 数 据 及 其 所 在 位 置 

终 采 用 KNN, 支持 向 量 机 (support vector machine，SVM)、 决 1.2 研究 方法 

策 树 (decision tree) 等 三 种 监督 方法 进行 分 类 ， 并 对 比分 析 三 本 文 首 先 对 无 人 机 获取 的 数据 进行 预 处 理 并 得 分 为 三 种 数 

种 监督 分 类 的 精度 结果 表明 数据 结合 的 提取 方法 能 够 获得 较 高 。 据 源 : a) 无 人 机 正 射 影像 直方 图 均衡 化 之 后 的 结果 作为 数据 源 

的 梯田 提取 结果 。 1; b) 用 DEM 数据 提取 相应 正 负 地 形 指数 PN (positive and 
"TM negative). 累计 曲率 AC (accumulative curvature ). JJ S£ S slope). 

! 实验 数据 和 研究 万 法 剖面 曲率 数 (profile curvature，PC)、 高 程 变异 系数 (coefficient 


1.1 实验 数据 of variation in elevation,CVE)、 粗 糙 度 (terrain roughness, TR), 
研究 区 为 中 国 甘肃 省 榆 中 地 区 龙泉 乡 黄土 丘陵 区 典型 旱 梯 ”山体 阴影 HS (hill shade) 等 7 种 地 形 指数 ， 并 对 7 种 地 形 指数 

田 ， 地 理 坐 标 范 围 东 经 104°10'58"~104°19'51", E Zi 进行 主 成 分 分 析 Cprincipal component analysis, PCA ) 降 维 处 理 ， 

35°34'4”~135°40'56”"。 数 据 于 2016 年 3 月 拍摄 ， 采 用 安 翔 动力 降 维 后 的 数据 作为 数据 源 2; c) 直方 图 均衡 化 后 的 无 人 机 影 

开发 的 AF1000 型 号 无 人 机 ， 搭 载 SONYA5100 相机 ， 单 幅 影 与 PCA 降 维 后 的 地 形 指数 结合 的 数据 作为 数据 源 3。 

像 面积 约 340mX500m, 影像 分 辩 率 0.05m。 拍摄 期 ， 风 力 小 于 其 次 基于 面向 对 象 方法 对 三 种 不 同 数据 源 进行 分 类 ， 

4 级， 天 气 上 晴朗 ， 能 见 度 高 ， 飞 行 采用 自动 起 飞 /规划 航线 飞行 。 对 分 类 精度 进行 评价 ， 进 行 面积 统计 。 研 究 的 技术 路 线 如 

/自动 降落 模式 ， 全 程 总 计 耗 时 约 24h。 数 据 处 理 采 用 Agisoft — 所 示 。 

photoscan 软件 ， 导 入 影像 ，POS 数据 及 控制 点 数据 ,将 整个 区 TR 可 | 
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Emp x 
ro Ul 


域 分 成 25 个 测 区 块 进行 处 理 ， 平 均 每 个 块 5000 张 左右 照片 ， E al BOE Lili 
Ll » oH | 1 | EARI 

对 每 个 块 逐一 进行 点 云 提取 和 立体 模型 建立 ,最 后 再 进行 全 并、 pp i cem oe] | | en | 
二 | 1 v | 

纹理 提取 ， 获 得 数字 表面 模型 (Digital Surface Model, DSM), 8 " | | [sse | o ona] | 

采用 CGCS2000 3 Gauss Kruger CM 105E 投影 。 由 于 研究 区 E E Sees | | | KETA Sup | 
L—ppewgGR ]—[ SIEERC]—[ rou ] | 0758 Y 

所 在 范围 在 该 段 时 间 几 乎 没有 植被 覆盖 ， 建 筑 用 地 较 少 ， 可 通 Ca 
| | 1 | Y | 

过 对 点 云 数据 进行 分 类 ， 分 为 地 面 点 和 非 地 面 点 两 类 ， 将 非 地 Imre: | 
| Ae] 二 x [8 | | 1 

面 点 作为 噪声 处 理 ， 通 过 滤波 后 得 到 吻 除 非 地 面 点 的 DEM 数 ZEE IEITICN | 
1 pe [m pm [8 [s HEC og pod | 

据 。DEM 分 辩 率 0.5m，30000 行 X24000 列 。 依 据 无 人 机 影像 IER IPS EE Lo 
bg a 

制作 的 1:500 地形 图 的 平面 和 高 程 精度 均 满 足 《1 : 5001: | tem 

1000 1: 2000 WE EIL" RE E PIA) UTE 1 : 500 平 图 2 本文 的 技术 路 线 图 

地 、 丘 陵 的 成 图 要 求 。 1.2.1 无 人 机 正 射 影像 直方 图 均衡 化 

本 文 截取 具有 不 同 特征 代表 性 的 三 块 1000pix*1000pix 数 高 分 辨 率 无 人 机 正 射影 像 地 物 景 观 的 结构 、 形 状 、 纹 理 和 


据 区 进行 试验 ， 如 图 1 所 示 。 区 域 一 的 梯田 区 有 部 分 山 丘 ， 建 ”细节 等 信息 较 突 出 ， 但 受 拍摄 天 气 状 况 影响 ， 易 出 现 对 比 度 不 
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足 ， 云 雾 模糊 等 现象 ， 为 特征 提取 增加 了 难度 。 通 常 采用 一 定 
的 方法 对 影像 进行 预 处 理 ， 以 增 大 对 比 度 ， 提 高 提取 精度 。 
直方 图 均衡 化 Chistogram equalization, HE) 03 是 常用 的 图 
像 对 比 度 增强 方法 。 HE 是 对 图 像 全 局 进行 调整 的 方法 , 在 某 些 
上 效果 较 差 : 自 适 应 直方 图 均值 化 Cadaptive histogram 
equalization, AHE) (4 算法 通过 计算 图 像 的 局 部 直方 图 ， 然 后 
重新 分 布 亮度 来 改变 图 像 对 比 度 ， 但 该 方法 在 部 分 区 域 有 过 度 


放大 图 像 中 相同 区 域 的 噪声 的 问题 ， 限 制 对 比 度 直方 图 均衡 化 


p 
Eu 


算法 Ccontrast limited adaptive histogram equalization, CLAHE) 
05 算 法 主要 是 通过 限制 AHE 算法 的 对 比 提高 程度 来 达到 的 ， 
对 每 个 小 区 域 都 必须 使 用 对 比 度 限 幅 ， 从 而 克服 了 AHE 的 过 
度 放大 噪声 的 问题 。 
在 本 文中 ， 结 合 CLAHE 方法 既 能 增强 影像 对 比 度 又 能 
服 噪声 过 度 放 大 的 优点 ， 对 无 人 机 正 射影 像 进 行 对 比 度 增强 和 
去 雾 预 处 理 ， 处 理 结果 作为 数据 源 1， 为 后 续 面 向 对 象 分 析 提 
供 良好 的 数据 支持 。 
1.2.2 地 形 指 数 的 计算 与 降 维 

梯田 的 田 面 识别 主要 是 对 田 坎 变 化 区 域 的 有 效 分 割 ， 目 视 
解 译 中 ， 作 为 地 形 指数 之 一 的 坡度 最 能 够 展现 该 变化 的 结果 。 
但 哪些 常用 的 地 形 指数 还 可 以 作为 特征 识别 的 参数 ， 需 要 进行 
界定 。 因 此 本 文 先 通过 ArcMap10.5 软件 ,通过 DEM 数据 计算 
出 相应 地 形 指数 : PN，S，AC，PC，CVE，TR，HS。 定 义 及 
计算 公式 如 下 : 

a)PN: 分 析 区 域内 所 有 栅 格 中 最 大 高 程 与 最 小 高 程 的 差 。 
WB, 为 局 部 最 大 高 程 值 ， 妃 ww 为 局 部 平均 高 程 值 。PN 的 计 
算 公式 如 下 : 


PN —H,,— H,, (1) 

b)S: 指 地 表单 元 陡 组 程度， 通常 把 披 面 的 生 直 高 度 HM 

水 平 距离 工 的 比 叫做 坡度 的 比 叫做 坡度 。 设 HS ERES S 的 
计算 公式 如 下 : 


S = arctan |=] (2) 

c)PC: 坡度 沿 地 面 最 大 坡 降 方向 的 地 面 高 程 变化 率 ， 本 文 
研究 区 的 坡度 以 度数 法 表示 。PC 的 计算 公式 如 下 : 

PC =Slop of S (3) 

d)AC: 地 面 曲率 是 对 地 形 表面 一 点 扭曲 变化 程度 的 定量 化 


度量 因子 ， 底 面 曲 率 在 垂直 和 水 平 两 个 方向 上 分 量 分 别称 为 剖 


率 (profile curvature) 和 平面 曲率 (planform curvature). 


累计 曲率 是 剖面 曲率 与 平面 曲率 的 差 值 。 设 Kp 为 剖面 曲率 ， 
Kc 为 平面 曲率 ， 则 AC 的 计算 公式 如 下 : 
AC = Kp— Kc (4) 
d)CVE: 反映 分 析 区 域内 地 表单 元 格 网 各 顶点 高 程 变化 的 
指标 ， 它 以 格 网 单元 顶点 的 标准 差 与 平均 高 程 的 比值 来 表示 ， 
设 SD (Standard Deviation) 为 标准 差 ，Z 为 领域 海拔 高 度 , CVE 
的 计算 公式 如 下 


CVE = SD/Z (5) 


FI 
RIR 


P LESS VE 信息 

e)HS: 山体 阴影 是 假想 一 个 光源 在 某 个 方向 和 某 个 太阳 高 

度 的 模拟 下 , 用 过 临近 像 元 的 计算 来 生成 一 副 0-255 的 灰 度 图 。 

假设 Hillshade < 0 ， 则 设 Hillshade 20 ，Zenith_rad 为 太阳 天 

顶 角 的 弧度 数 , Slope. rad 为 菜 点 的 坡度 弧度 数 , Azimuth _ rad 

为 太阳 光线 方向 角 的 弧度 数 ， spect_rad 为 某 一 点 的 坡 向 弧度 
数 ， 则 HS 的 计算 公式 如 下 : 


cos Zenith rad * 
cos Slope rad 


HS = 255.0 * '|sin Zenith rad * (6) 
sin Slope rad * 
: pan _rad i 
Aspect rad 


ATR: 地 表单 元 的 曲面 面积 S 与 其 在 水 平面 上 的 投影 面 


之 比 。 计 算 公式 如 下 : 


TR = 1/cos S (7) 

在 包含 高 程 、 坡 度 等 的 多 波段 栅 格 中 ， 地 形 指数 通常 由 高 
程 数据 直接 或 间接 得 出 ， 数 据 元 余 现象 较为 明显 。 为 了 减少 在 
面向 对 象 分 类 时 的 数据 元 余 , 本 文采 用 主 成 份 分 析 法 09PCA 进 
行 数据 降 维 后 ， 作 为 数据 源 2。 
1.2.3 面向 对 象 分 类 及 评价 
本 文 利用 面向 对 象 方法 对 梯田 的 田 面 进 行 提 取 。 首 先 对 三 
类 数据 源 进行 多 尺度 分 割 的 最 优 尺度 选择 ， 将 选择 的 结果 作为 
专题 数据 ， 在 eCognition Developer 9.0 软件 中 进行 梯田 的 田 面 
提取 .面向 对 象 分 类 采用 KNN,SVM, 决策 树 等 三 种 监督 方法 ， 
对 比分 析 三 种 监督 分 类 的 精度 。 
a) 最 优 尺度 分 割 。 对 于 基于 面向 对 象 分 析 的 遥感 图 像 地 物 
提取 而 言 ， 图 像 分 割 是 第 一 步 ， 图 像 分 割 的 质量 直接 决定 着 后 
续 的 面向 对 象 分 类 提取 梯田 的 精度 073。 近 年 来 ， 图 像 分 割 方法 
的 研究 主要 概括 为 基于 图 论 ， 聚 类 ， 分 类 等 算法 08 29。 目 前 基 
于 超 像素 的 图 像 分 割 法 较为 流行 202， 如 SLIC (Simple Linear 
Iterative Clustering )??!, SEEDS Superpixels Extracted via Energy- 


Driven Sampling) P4, LSC (Linear Spectral Clustering ) 251, Mean 


Shift?9!, Marker-based WatershedP27、Graph-CutsP23 等 等 。 基 于 
超 像 素 法 对 影像 进行 分 制 ， 能 有 效 利 用 影像 的 光谱 、 纹 理 、 形 
状 等 特征 ， 能 获得 较为 均 质 的 对 象 。 但 是 ， 目 前 的 超 像素 分 割 
不 能 产生 细节 保留 和 过 分 割 结 果 。 大 多 数 情况 下 ， 超 像素 分 割 
产生 大 量 的 过 分 割 碎片 ， 造 成 大 量 数 据 元 余 ， 在 图 像 分 类 和 

标 检 测 中 ， 元 余 的 分 割 碎片 很 难 提 高 效率 和 重 棒 性 ， 而 且 还 增 


加 计算 负载 09。 因 此 , 需要 在 分 割 的 基础 上 进行 进一步 的 处 理 。 
分 割 后 处 理 阶 段 通常 是 根据 分 割 目的 ， 对 分 割 后 的 结果 再 次 进 
行 处 理 ， 例 如 区 域 合 并 及 区 域内 滤波 去 噪 等 。 

前 较为 流行 的 FNEAP9 多 尺度 分 割 就 是 基于 超 像素 初始 
分 割 后 ， 再 进行 区 域 合并 获取 分 割 对 象 。 多 斥 度 分 割 法 是 一 种 
自 下 而 上 的 分 割 方法 ， 从 单个 像 元 开始 ， 向 上 逐渐 合并 成 较 大 
的 对 象 ， 直 到 满足 所 设置 的 分 割 参 数 条 件 。 与 基于 单 像 元 光谱 
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特性 的 影 


像 信息 提取 方法 相 比 ， 多 


尺度 分 割 能 够 更 充分 地 利用 


高 分 辩 率 遥感 影像 的 信息 。eCognition Developer9.0 软件 提供 


FNEA 多 


解 译 人 员 和 


R 


据 研究 


限制 了 高 分 辩 率 遥感 影像 的 


趋势 B0。 


mI Zh ARE RI 


究 方 法 


JÉ 


E3 
据 最 小 


将 每 个 


aN 
as 


ES 


[5 


尺度 分 割 算法 ， 所 涉及 主要 参数 
致 度 3 个 ， 尤 其 是 与 分 割 对 象 大 小 密切 机 
区 域 、 目 标 并 结合 以 往 经 验 设 定 。 这 极 大 地 
自动 解 译 。 借 鉴 人 工 试 错 方法 对 分 
已 成 为 当前 的 一 种 新 


价 ， 实 现 分 割 参 数 自动 选择 


自 适应 地 选取 较为 合适 的 分 割 
度 分 割 技术 中 的 关键 问 
法 ， 在 尺度 集 模型 中 进行 全 
督 尺度 集约 简 , 得 到 最 优 分 
用 层次 迭代 优化 的 区 域 合并 方法 ， 构 建 了 区 域 层 次 结构 ， 
| 影像 的 多 尺度 表达 ， 即 
区 域 合并 过 程 完整 的 记录 下 来 ， 并 记录 整个 过 程 中 产生 
4 层次 关系 ， 
集 模型 可 以 反 算 
度 参数 的 问题 。 
风险 贝 叶 


题 


IRE. 


斯 决策 规则 进行 基于 全 
简 ， 最 后 再 基于 局 部 演化 分 析 进 行 


。 本 文 应 | 


区 域 进行 
UE RC 
区 域 合并 的 同时 进行 了 


尺度 是 
JHH 


包括 
目 关 的 


高 分 辨 率 遥 感 


er^ A 
忠文 等 [ 


RE HAR 


的 
紧 


hinaXi 


C 
BALA, LAMB UU AS e M i CUR 


练 集中 与 之 最 为 相似 的 前 天 个 数据 , 则 该 测试 数据 对 应 的 类 别 
就 是 个 数据 中 出 现 次 数 最 多 的 那个 分 类 。 计算 对 象 间 不 同 特 
征 值 之 闻 的 距离 作为 各 个 对 象 之 间 的 非 相 似 性 指标 ， 避 免 了 对 


尺度 参数 ， 需 


区 


"0 的 最 


化 分 析 并 进行 非 


该 研究 结合 FNEAD9 算 法 ， 


尺度 集 模型 


。 尺 度 集 模 型 可 


MER 


| 结果， 从 而 
全 局 演化 分 析 ， 
局 演化 分 析 的 尺度 集约 
尺度 集约 简 。 选 择 欠 分 割 的 


引 。 


此 通过 尺 


要 


像 
优 


TVs 
m. 


以 


度 


解决 调整 分 


代价 权重 C-1 时 ， 尺 度 集 达到 最 佳 约 简 ， 选 择 其 为 最 佳 分 割 


度 。 


bo 分 类 特征 筛选 。 多 尺度 分 割 完 成 之 后 ， 本 文 先 基于 
田 和 非 梯田 样本 对 象 ， 
FE 计算 ， 计 算 方 法 如 表 1。 对 
图 统计 直观 分 析 ， 筛 选 可 分 性 较 强 的 


一 的 三 种 数据 源 ， 
对 样本 对 象 的 各 波段 进行 分 类 特 生 
得 到 的 分 类 特征 进行 箱 线 


Hi 


LAEE H 


O 


区 


定量 的 梯 


38] 


依 


R 


域 


特征 作为 分 类 特征 。 
表 1 AGENS 
特征 说 明 
平均 值 Mean 影像 对 象 的 所 有 像素 的 图 层 值 的 平均 值 
标准 差 StdDev 影像 对 象 的 所 有 像素 的 图 层 值 计 算得 到 标准 差 
m f 影像 对 象 的 图 层 数量 除 以 包含 光谱 信息 的 图 层 
亮度 Brightness 
平均 值 的 总 和 
长 宽 比 
协 方差 窍 阵 的 特征 值 的 比值 ， 较 大 的 作为 分 子 
Length/Width 
形状 指数 影像 对 象 的 边界 长 度 除 上 它 的 面积 的 平方 根 的 


Shape Index 


AH 


n 


ÆI KE 


GLCM 


4 fi 


相关 性 CCorrelation) 
对 比 度 (Contrast) 


Ii (Entropy) 


均 质 性 (Homogeneity) 


c) 监 督 分 类 。 应 用 分 类 特 生 


SVM， 决 策 树 监督 分 类 。 


KNN 算法 


EX 


RUE Ap S BE Sod 


fT KNN, 


是 一 种 非 参数 的 ,基于 统计 的 分 类 算法 。 主 体 思 


想 : 在 训练 集中 数据 和 标签 已 知 的 情况 下 ， 输 入 测试 数据 ， 将 


测试 数据 的 特征 与 训练 集中 对 应 的 特 生 


进行 相互 比较 ， 找 到 训 


象 之 间 的 匹配 问题 , 距离 计算 方法 一 般 使 用 欧 氏 距离 ( 式 (8))。 
对 于 本 文 的 关于 梯田 与 非 梯田 的 二 分 类 问题 ，K 值 应 选择 为 奇 


dose Bre cs 
ka 


SVM 算法 ， 最 早 由 Cortes 和 Vapnik 提出 B1， 是 一 个 凸 二 


(8) 


次 规划 求解 问题 ， 即 对 线性 不 可 分 的 两 类 问题 ， 其 最 优 的 分 类 
形式 为 : 把 学 习 样本 非 线性 映射 到 高 维特 征 空间 ， 在 高 维特 征 
空间 构造 线性 判别 函数 ， 将 原样 本 空间 中 非 线性 可 分 的 问题 转 


化 为 在 高 维特 征 空间 线性 可 分 的 问题 ,使 得 两 类 无 错误 地 分 开 ， 
使 两 类 的 分 类 间隙 最 大 。 在 小 样本 、 非 线性 情况 下 ， 具 有 和 较 
好 的 泛 化 性 能 。 本 文 基于 径 向 核 函 数 (radial basis function,RBF )， 
遗传 算法 (genetic algorithms, GA) 寻找 最 优 参数 〈 惩 罚 


对 象 


采用 
系数 c 与 核 函 数 半径 g)， 对 多 尺度 分 割 后 的 梯田 与 非 梯 | 
建立 分 类 模型 。 

决策 树 是 一 种 基本 的 分 类 与 回归 方法 。 决 策 树 模型 呈 树 
结构 ， 在 分 类 问题 中 ， 表 示 基 于 特征 对 实例 进行 分 类 的 过 程 。 
决策 树 学 习 算 法 包含 特征 选择 、 决 策 树 的 生成 与 前 枝 过 程 。 分 
类 与 回归 树 (classification and regression tree, CART) 算法 B3 选 
择 特征 采用 基尼 指数 (Gini coefficient)。 将 训练 数据 集 ( 训 练 
样本 ) 划分 为 测试 变量 和 目标 变量 ， 通 过 对 这 两 变量 的 循环 分 
析 形 成 二 叉 决 策 树 。 基 于 对 象 的 决策 树 分 类 在 地 物 提 取得 到 广 
泛 应 用 B3 约 ， 能 更 好 的 挖掘 图 像 信 息 ， 是 图 像 分 类 从 一 种 纯粹 
的 认识 过 程 转化 为 一 种 县 有 先 验 知 识 的 再 认识 过 程 B9 。 

d) 分 类 精度 评价 。 精 度 评 价 是 指 比较 实地 数据 与 分 类 结果 ， 
以 确定 分 类 过 程 的 准确 程度 。 本 文通 过 目 视 解 译 获得 分 类 评价 
样本 。 在 arcmap10.5 里 绘制 接近 于 实际 梯田 的 真 值 图 ， 并 且 统 
计 了 田 面 的 实际 面积 。 基 于 目 视 解 译 获 得 的 评价 样本 ， 采 用 
eCognition Developer9.0 软件 中 提供 的 “Error Matrix based on 
TTA MASK” 对 分 类 结果 进行 评价 。 最 常用 的 精度 评价 方法 是 
混淆 和 矩阵 B33 (Error Matrix) 方法 。 设 样本 总 类 别 数 为 k ， 样 
本 总 数 为 Y ,测试 样本 中 i 类 的 样本 总 数 和 N,, 计算 公式 如 式 (9) 
所 示 。 实 际 分 类 中 被 分 为 j 类 的 样本 总 数 N,; 的 计算 公式 如 式 

(100 所 示 。 


NS 


(9) 


(10) 


其 中 : Ns 为 测试 样本 中 本 应 该 分 为 i 类 而 实际 分 类 结果 中 却 被 
分 为 j 类 的 样本 数目 。 从 混淆 矩阵 可 以 计算 出 各 种 精度 统计 值 ， 
如 生产 者 精度 (producer’s accuracy, PA) P581, fj PTS EE (user’s 
accuracy，UA)， 总 体 精 度 (overall accuracy; OA), Kappa 系 
数 户 ! 等 ， 计 算 公式 分 别 如 式 〈11) ~ (140. 所 示 。 
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NYN, -$o *N,;) 


K= 


k 
N? -Y (N, *N,,) 
i=1 


2 ”结果 与 讨论 


2.1 CLAHE 处 理 结果 


本 文 对 无 人 机 正 射 影像 采用 CLAHE 算法 进行 增强 (图 


an 


(12) 


(13) 


(14) 


3), 


对 比 增强 前 后 可 以 得 出 ,增强 后 的 影像 数据 | 
为 进一步 研究 提供 良好 的 数据 源 。 
2.22 地形 指数 计算 及 降 维 结果 


坎 轮廓 更 加 清晰 ， 


Chi 
无 人 机 影像 nul 指数 台 


SX bata 


BALA, F: a 
表 4 区 域 3 地 形 指数 相关 性 分 析 
波段 PN AC CVE S PC TR HS 
PN 1.000 
AC 0.031 1.000 
CVE 0.978 -0.010 1.000 
S 0.894 -0.003 0.903 1.000 
PC 0.229 0.002 0.234 0.268 1.000 
TR -0.000 0.000 -0.000 -0.001 -0.001 1.000 
HS -0.177 0.003 -0.180 -0.135 -0.052 1.000 


利用 


进行 


ArcMap10.5 软件 求 出 三 
PCA 分 析 ， 得 出 相关 系数 矩阵 ， 


块 样 区 的 七 种 地 形 指 数 ， 并 


如 表 2，3，4 所 示 。 观 察 


三 个 相关 性 分 析 表 ,得 出 ,PN 与 CVE 的 相关 系数 分 别 为 0.980， 
0.975, 0.978. PN 5E S 的 相关 系数 分 别 为 0.924, 0.910, 0.894。 
CVE 5 S 的 相关 系数 分 别 为 0.939，0.929，0.903。 设 相关 系数 


闷 值 为 0.9， 可 以 确定 CVE, PN, S 
ANR. X7 种 地 形 指数 进行 PCA 主 成 份 分 析 , 得 


之 间 具 有 强 关 联 性 ， 信 息 
出 5 种 相关 


性 弱 的 数据 组 合 。 实 验 结果 为 S，PC，AC，HS，TR 这 5 种 地 


形 指数 为 主要 成 份 ， 如 图 4 所 示 。 其 中 S1，PC1，AC1，HS1， 

TRI 代表 区 域 一 的 ， 依 此 类 推 。 
表 2 区 域 1 地 形 指数 相关 性 分 析 

波段 ”PN AC CVE S PC TR HS 

PN 1.000 

AC 0.05 1.000 

CVE 0.980 -0.010 1.000 

S 0.924 -0.007 0.939 1.000 

PC 0.71 0.004 0.175 0.191 1.000 

TR -0.001 -0.000 -0.001 -0.001 0.001 1.000 

HS  -0.328 0.005 -0.330 -0.332 -0.040 -0.000 1.000 
表 3 区 域 2 地 形 指数 相关 性 分 析 

波段 ” PN AC CVE S PC TR HS 

PN 1.000 

AC 0.034 1.000 

CVE 0.975 -0.009 1.000 

S 0.910 -0.003 0.929 1.000 

PC 0237 -0.001 0241 0256 1.000 

TR -0.000 -0.001 -0.000 -0.000 0.000 1.000 

HS -0.016 0.005 -0.003 0.000 0.030 -0.002 1.000 


(a 


v 


区 域 1 原始 图 


(c) 区 域 2 原始 图 


区 域 1 


-0.000 


CLAHE 后 


(e 


— 


区 域 3 原始 图 (f) 
图 3 CLAHE 前 后 对 比 图 
23 分 类 及 结果 评价 
2.3.1 多 尺度 分 割 

ETE 
将 形状 因子 为 0.5， 紧 致 度 因 子 为 0.8， 
出 最 佳 分 割 尺 度 ， 其 次 
Developer9.0 中 ， 作 为 专题 数据 ， 


究 区 域 的 三 类 数据 源 ， 首 先 应 用 胡 上 忠文 等 I%11 方 法 ， 
建立 相应 的 尺度 自 
将 分 析 结 果 导 入 到 eCognition 
进行 多 尺度 分 割 , 如 图 5 所 


区 域 3 CLAHE 后 


找 


Yit 


示 。 对 比 基 于 三 种 数据 源 的 分 割 结果 ， 分 析 
d), 
线 区 域 。 


区 域 2 图 Sf，g，h) ,区 域 3 RI 5j, k, D 的 红 ， 蓝 


a) 基于 数据 源 1 (图 5b，f，j) 的 分 割 。 


区 域 1 图 5b, c, 


Ani » 


红 框 和 蓝 框 区 域 


边缘 光谱 差异 较 小 ， 履 盖 


都 出 现 欠 分 割 。 欠 分 割 的 原因 : Bil 
Y, HRED EKI), I 
导致 无 法 分 割 。 


坎 过 窗 〈 蓝 框 区 域 ) 等 的 干扰 ， 
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p à; II Ga qp tec ^ e. d 
5 ON Aca HS3 TRÀ 
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图 5 研究 区 域 1-3 的 三 类 数据 源 的 多 尺度 分 割 最 优 效果 图 


到 

T as 
-E 
Am 


es vtr AA )) 
(i) 区 域 三 的 真 值 


b) 基 于 数据 源 2 (图 Sc. g. O 的 分 割 。 在 区 域 一 、 三 《图 分割 处 于 欠 分 割 状 态 。 而 且 区 域 2《〈 图 5g) 分 割 结果 较 差 ， 边 
Sc, k) 的 红 框 和 蓝 框 区 域 分 割 完好 ， 边缘 清 晰 可 见 。 与 基于 数 。” 缘 模 糊 难 辨 。 分 割 效果 差异 原因 : 在 光谱 差异 不 明显 的 梯田 边 
据 源 1 区 域 1、 区域 3( 图 Sb、j) 分 割 比较 〈 红 框 和 蓝 框 区 域 )， 缘 区 域 ， 地 形 指数 信息 却 有 很 好 的 可 分 性 ， 能 良好 地 为 梯田 分 
都 能 很 好 地 区 分 出 梯田 边缘 ， 但 在 绿 框 区 域 ， 基 于 数据 源 2 的 ”人 割 提 供 数据 支持 ， 但 是 从 区 域 2 的 地 形 指数 (图 4、S2、PC2、 
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Chi 
RAM BAR, $: aiant aU 结 从 的 


AC2、HS2、TR2) 分 析 来 看 ， 区 域 2 (图 5g) 的 梯 


于 田 面 ”特征 。 随 机 选取 梯田 与 非 梯田 样本 ， 尽 量 使 样本 保持 平衡 。 对 


TE, AEDE EMRE) 等 原因 ， 地 形 指数 已 经 不 能 充 “ 于 基于 数据 源 1， 数 据 源 2 的 分 类 ， Sm TRUE Hen. 
分 表达 梯田 信息 。 也 可 对 比 绿 框 (图 5, j O 区 域 分 析 , 发 现 ” 算 ( 表 2)， 然 后 参与 分 类 。 但 是 对 数据 源 3， 由 于 特征 数量 过 
对 于 平坦 的 田 面 ， 当 好 几 块 田 处 于 同一 个 平面 ， 虽 然 光谱 差异 多、 数据 元 余 ， 需 要 进行 特征 优化 。 统 计 样 本 的 特征 ， 筛 选 出 


大 ， 但 地 形 指数 的 可 区 分 行 就 相对 较 弱 。 能 区 别 梯田 与 非 梯 田 较为 明显 的 特征 信息 , 制作 箱 线 图 (图 6)。 
c) 基于 数据 源 3( 图 5d、h、1) 的 分 割 。 在 区 域 1~3 的 红 “结果 表明 , R 波段 的 均值 (Mean R), S 波段 的 均值 (Mean S), 

框 和 复 框 中 ， 分 割 结 果 边 缘 清 晰 可 见 ， 将 基于 数据 源 1 和 数据 G, B, TR, HS 波段 的 标准 差 (Standard deviation G, B, TR, 

源 2 的 优势 结合 ,充分 地 得 到 相对 较为 接近 真 值 图 的 分 割 结果 。 HS), S 波段 的 纹理 特征 CGLCM correlations S, GLCM 

2.3.2 监督 分 类 homogeneity S, GLCM deviation S) 区 分 度 较 大 ,适合 作为 分 类 
1) 分 类 特征 提取 ”对 多 尺度 分 割 后 的 对 象 进行 特征 提取 ， 特征 参数 。 

分 别 求 对 象 的 R，G，B，S，PC，AC，HS，TR 等 各 个 波段 的 
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图 6 ”多 尺度 分 割 对 象 的 特征 信息 
2) 分 类 与 评价 ”本 文 首先 基于 区 域 一 的 三 种 数据 源 进 行 (c) 基 于 数据 源 3 的 三 种 监督 分 类 ， 都 不 同 程度 地 出 现 错 
KNN, SVM, CART 的 监督 分 类 研究 ， 其 中 图 7 是 区 域 一 的 真 。” 分 现象 ， 其 中 黑 框 区 域 都 出 现 错 分 现象 ， 图 10d 内 紫 框 区 域 ， 


ER, 图 8~10 为 分 类 结果 , 红色 区 域 为 梯田 的 田 面 , 黄色 区 域 E 10d 蓝 框 区 域 出 现 错 分 。 错 分 的 原因 : 分 类 特征 的 选取 和 训 
为 非 梯 田 区 域 。 分 类 精度 评价 如 表 5 所 示 。 练 样本 的 选择 差异 及 分 类 方法 所 致 。 
a) 分 类 结果 定性 分 析 如 下 : 总 结 以 上 分 析 结 果 发 现 ， 影 响 分 类 结果 的 首要 因素 是 多 尺 


(a) 基于 数据 源 1 的 三 种 监督 分 类 , 在 蓝 ， 黑 ， 紫 框 内 区 ” 度 分 割 的 效果 ， 基 于 数据 源 1 的 多 尺度 分 割 ， 在 梯田 边缘 光谱 
域 ， 都 出 现 错 分 现象 。 错 分 原因 : 蓝 ， 紫 框 内 区 域 ， 在 多 尺度 。 ”变化 不 明显 的 区 域 会 容易 出 现 欠 分 割 ,从 而 导致 分 类 结果 出 错 ; 
分 割 时 ， 因 光谱 特征 与 周边 对 象 差异 不 大 而 出 现 从 分 割 ， 而 黑 ”基于 数据 源 2 的 多 尺度 分 割 ， 在 光谱 差异 很 大 ， 但 地 形 指数 差 
框 区 域 可 能 因 样 本 特征 的 选取 不 完整 或 分 类 方法 所 致 ， 异 不 明显 的 地 方 出 现 欠 分 割 ， 从 而 导致 分 类 结果 出 错 ; 基于 数 

(b) 基于 数据 源 2 的 三 种 监督 分 类 , 在 黑 ， 紫 ， 绿 框 内 区 ” 1808 3 的 多 尺度 分 割 ， 参 考 图 5d、h、1， 多 尺度 分 割 结果 都 优 
域 出 现 错 分 现象 。 错 分 原因 : 由 于 地 形 指数 差异 不 明显 而 导致 于 基于 数据 源 1、2， 但 分 类 结果 出 现 不 同 程度 的 差异 ， 是 因为 
多 尺度 分 割 处 于 欠 分 割 状 态 。 而 其 他 错 分 的 原因 是 地 形 指数 特 ”监督 分 类 的 影响 。 至 此 ， 本 文 相对 充分 地 验证 了 ， 基 于 数据 源 
征 并 没 充分 表达 梯田 特 和 3 作为 梯田 的 田 面 提取 的 数据 源 最 为 合适 。 对 于 基于 三 种 监督 
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分 类 的 影响 ， 本 文 基于 分 类 精 
b) 分 类 结果 定量 分 析 如 下 : 
对 比分 析 表 5 中 每 一 种 分 类 方法 


度 评价 。 


度 评价 所 


ÆJ 


产 者 精度 和 
类 方法 ， 基 于 数据 源 1 分 类 


] 户 精度 角 


VE 


度 分 
取 的 生产 者 精 


日 


析 。 


ck p 


如 


基于 三 种 数据 源 的 分 类 精 


MH KNN 分 
青 度 ， 用 户 精 度 分 别 


对 象 分 


CART 的 三 种 监督 分 类 结果 如 图 11 所 示 。 分 类 精度 记 
所 示 。 对 


分 类 法 ， 


系数 为 0.872，0.732。 


Chi 
薛 牡 丹 ， 等 : 无 人 机 影像 s IDEE 


Wr. 多 尺度 分 


比分 析 表 5 


割 结果 如 


NISNI 


图 5h, 1 所 示 。 基 于 KNN、SVM、 


EF 价 如 表 6 


区 域 1 的 


< 


于 数据 源 3 的 精度 评价 ，KNN 


总 精度 和 Kappa 系数 为 0.887，0.739; SVM 分 类 法 总 精 
度 和 Kappa 系数 为 0.887, 0.765; CART 分 类 法 总 精度 和 Kappa 


显然 SVM 分 类 法 优 于 KNN 法 和 CART 


为 0.759，0.908 。 基 于 数据 源 2 分 类 提取 的 生产 者 精度 ， 用 户 ” 法 。 同 理 对 比分 析 表 8， 区 域 2、3 的 基于 数据 源 3 的 三 种 分 类 
精度 分 别 为 0.744，0.771 。 基 于 数据 源 3 分 类 提取 的 生产 者 精 。 评价 。 对 于 区 域 2，SVM 分 类 法 优 于 KNN 法 和 CART 法 。 对 
度 ， 用 户 精度 分 别 为 0.776，0.908。 显 然 基于 数据 源 3 的 分 类 ”于 区 域 3，CARIT 分 类 法 优 于 KNN 法 和 SVM 法 。 
提取 梯田 的 错 分 最 少 。 同 理 分 析 应 用 SVM，CART 分 类 方法 ， 从 实验 整体 结果 分 析 , 基于 SVM 分 类 法 优 于 基于 KNN 法 
基于 数据 源 3 分 类 提取 的 生产 者 精度 ， 用 户 精 度 几乎 都 高 于 基 M CART 法 。 从 算法 原理 的 角度 分 析 ，KNN 法 简单 ， 有 效 , 重 
于 数据 源 1，2 新 训练 的 代价 较 低 , 但 计算 量 较 大 ,CART 法 易于 理解 和 解释 ， 
(b) 从 总 体 精度 和 Kappa. 系数 角度 分 析 。 如 应 用 KNN 的 。” 能够 同时 处 理 数据 型 和 常规 型 属性 。 但 对 于 各 类 别 样本 数量 不 
分 类 方法 ， 基 于 数据 源 1 提取 为 梯田 的 总 体 精度 ，Kappa 系数 一致 的 数据 ， 在 CART 当中 ,信息 增 益 的 结果 偏向 于 有 具有 更 多 
分 别 为 0.875，0.730。 基 于 数据 源 2 提取 为 梯田 的 总 体 精度 数值 的 特征 。SVM 法 可 以 解决 小 样本 情况 下 分 类 ,可 以 提高 泛 
Kappa 系数 分 别 为 0.809，0.600。 基 于 数据 源 3 提取 为 梯田 的 ”化 性 能 ， 解 决 高 维 问题 ， 非 线性 问题 。 结 合 本 文 的 数据 特点 ， 
总 体 精度 ，Kappa 系数 分 别 为 0.887，0.739。 显 然 基 于 数据 源 3 ”研究 区 域 的 梯田 与 非 梯田 样本 数量 有 限 ， 样 本 特征 多 样 ， 较 为 
的 提取 梯田 的 分 类 精度 高 于 基于 数据 源 1、2。 同 理 分 析 应 用 ”适合 基于 SVM 法 进行 分 类 提取 。 本 文 最 终 以 SVM 的 分 类 结 
SVM, CART 分 类 方法 ， 基 于 数据 源 3 的 总 体 精度 ，Kappa 系 果 ， 对 梯田 的 面积 进行 了 统计 ， 如 表 7 所 示 。 与 目 视 解 译 获得 
数 都 高 于 基于 数据 源 1、2。 因 此 ， 将 正 射 影像 数据 与 地 形 指 数 ”的 面积 差距 相对 较 大 。 其 主要 原因 是 基于 SVM 法 的 核 函 数 及 
相 结 合 的 梯田 的 田 面 提取 方案 优 于 单一 数据 源 提取 方案 。 参数 设置 的 不 够 合理 ， 同 时 目 视 解 译 存在 一 定 的 不 可 避免 的 误 
根据 区 域 1 的 研究 结果 , 对 区 域 2、3 的 数据 源 3 进行 面向 — 差 也 影响 着 参考 数据 。 


bo 1 的 真 值 图 


基于 数据 源 1 分 类 


基于 数据 源 1 分 类 


(b) 


g 


F 数据 源 1 分 类 


8 区 域 1 的 三 种 数据 源 的 KNN 分 类 对 比 


区 | 


(c) 基于 


数据 


2 分 类 
9 区 域 1 的 三 种 数据 源 的 SVM 分 类 对 比 


(c) 4&7 


数据 源 2 分 类 
10 区 域 1 的 三 种 数据 源 的 CART 分 类 对 比 


源 2 分 类 


(d) 基 


(Cd) 基 


(d) 基于 数据 源 3 分 类 


于 数据 源 3 分 类 


于 数据 源 3 分 类 
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(a) 区 域 2 真 值 


(e) 区 域 3 真 值 图 (DKNN 分 类 


图 11 


区 域 2.3 基于 数据 源 3 的 监督 分 类 对 比 图 


表 5 区 域 1 的 分 类 精度 评价 


图 (b)KNN 分 类 


合作 期 刊 
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(c)SVM 分 类 (d)CART 分 类 


(g)SVM 分 类 (h)CART 分 类 


表 7 区 域 1~3 的 面积 统计 


分 类 方法 及 精度 评价 
分 类 方法 ”精度 评价 


数据 源 1 ”数据 源 2 


数据 源 3 


工 NT T NT T NT 


生产 者 精度 0.759 0.95 0.744 0.852 0.776 0.944 
用 户 精度 ^ 0.908 0.859 0.771 0.838 0.908 0.868 
TM 总 精度 0.875 0.809 0887 
Kappa 系数 0.730 0.600 0.739 
生产 者 精度 0.757 0.935 0.807 0.857 0.876 0.894 
户 精 度 — 0.883 0.894 0.792 0.866 0.848 0.919 
Pos 总 精度 0.865 0.838 0.887 
Kappa 系数 0.710 0.662 0.765 
生产 者 精度 0.732 0.954 0.818 0.840 0818 0.908 
户 精 度 ^ 0.912 0.854 0.774 0.879 0.858 0.886 
id 总 精度 0.867 0.832 0.872 
Kappa 系数 0.711 0.625 0.732 
表 6 区 域 2、3 的 基于 数据 源 3 的 分 类 评价 
分 类 方法 及 精度 评价 区 域 二 区 域 三 
分 类 方法 精度 评价 T NT T NT 
生产 精度 0.870 0.770 0.702 0.834 
户 精度 0.857 0.787 0.853 0.768 
S 总 精度 0.829 0.811 
Kappa 系数 0.644 0.621 
生产 精度 0.897 0.778 0.761 0.874 
户 精 度 0.854 0.839 0.880 0.807 
Md 总 精度 0.848 0.812 
Kappa 系数 0.682 0.626 
生产 精度 0.770 0.856 0.990 0.775 
户 精 度 0.885 0.720 0.832 0.885 
Mi 总 精度 0.805 0.855 
Kappa 系数 0.608 0.704 


区 域 视 解 译 面 积 SVM 分 类 面积 ES 
区 域 一 98484.50 93811.50 -4673.00 
区 域 二 147601.00 155128.75 +7527.75 
区 域 三 137070.50 139780.25 +2709.75 
3 ARE 


源 的 结合 的 面向 对 象 梯田 信息 提取 。 研 究 表明 基于 无 人 机 正 射 
影像 与 地 形 指 数 结合 的 数据 ， 充 分 利用 了 影像 的 光谱 、 纹 理 、 
形状 信息 、 以 及 地 形 指数 信息 ， 结 合 SVM 分 类 方法 的 面向 对 
象 提取 梯田 ， 精 度 达 到 较 好 的 效果 。 


本 文 基于 高 分 辨 率 无 人 机 正 射 影像 、 地 形 指数 和 两 种 数据 


T 


本 文 还 有 许多 需要 改进 之 处 : a) 在 基于 SVM 分 类 方法 的 面 


Ed 


向 对 象 梯田 提取 中 ， 参 数 择优 应 用 GA 算法 ， 但 是 否 存在 更 合 


T 


BRIR, AREP RR bE 
了 箱 线 图 进行 了 统计 分 析 ， 并 未 探索 更 为 合理 的 优化 方法 ;c) 能 


本 特征 优化 只 应 用 


和 否 将 多 种 分 类 方法 结合 起 来 使 用 ， 即 集成 机 器 学 习 方 法 ， 或 者 


应 | 


深度 学 习 等 方法 进行 分 析 ， 有 待 进行 进一步 探索 。 
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